Fast rotating condensates in an asymmetric harmonic trap 
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We investigate the effect of the anisotropy of a harmonic trap on the behaviour of a fast rotating 
Bose-Einstein condensate. Fast rotation is reached when the rotational velocity is close to the 
smallest trapping frequency, thereby deconfining the condensate in the corresponding direction. We 
characterize a regime of velocity and small anisotropy where the behaviour is similar to the isotropic 
case: a triangular Abrikosov lattice of vortices, with an inverted parabola profile. Nevertheless, at 
sufficiently large velocity, we find that the ground state does not display vortices in the bulk. We 
show that the coarse grained atomic density behaves like an inverted parabola with large radius 
in the deconfined direction, and keeps a fixed profile given by a Gaussian in the other direction. 
The description is made within the lowest Landau level set of states, but using distorted complex 
coordinates. 

PACS numbers: 



Vortices appear in many quantum systems such as su- 
perconductors and superfluid liquid helium. Rotating 
atomic gaseous Bose-Einstein condensates constitute a 
novel many body system where vortices have been ob- 
served [l[ and various aspects of macroscopic quantum 
physics can be studied. In a harmonically trapped con- 
densate rotating at a frequency close to the trap fre- 
quency, interesting features have emerged, presenting a 
strong analogy with quantum Hall physics. In the mean 
field regime, vortices form a triangular Abrikosov lattice 
and the coarse grained density approaches an inverted 
parabola [!, H, 0] • At very fast rotation, when the num- 
ber of vortices becomes close to the number of atoms, 
the states are strongly correlated and the vortex lattice 
is expected to melt [3]. In the mean field regime, Ho 
[H observed that the low lying states in a symmetric 2D 
trap are analogous to those in the lowest Landau level 
(LLL) for a charged particle in a uniform magnetic field. 
This analogy allows a simplified description of the gas by 
the location of vortices: the wave function describing the 
condensate is a Gaussian multiplied by an analytic func- 
tion of the complex variable z = x + iy. The zeroes of 
the analytic function are the location of the vortices. It 
is the distortion of the vortex lattice on the edges of the 
condensate which allows to create an inverted parabola 
profile @, i, 0, i for the coarse grained atomic density 
in the LLL. 

The experimental achievement of rotating BEC in- 
volves anisotropic traps. An anisotropy of the trap can 
drastically change the picture in the fast rotation regime. 
In this case, the condensate becomes very elongated in 
one direction and forms a novel quantum fluid in a nar- 
row channel. The investigation of the vortex pattern has 



been performed for an infinite strip which corresponds to 
the situation where the rotational frequency has reached 
the smallest trapping frequency [ToL [Txj] , and for an elon- 
gated condensate (l2l H3 . [l4j . As pointed out by Fetter 
[Tj| , the description of the condensate can still be made 
in the framework of the lowest Landau level, defined by 
an anisotropic Gaussian, multiplied by an analytic func- 
tion of x + i/3y, where (3 is related to the anisotropy of 
the trap and the rotational frequency. We are going to 
characterize a regime of fast rotation where there are 
no vortices in the bulk of the condensate and show that 
the coarse grained density profile is very different from 
the isotropic case: the behaviour is an inverted parabola 
with large expansion in the deconfined direction, while 
the extension remains fixed in the other direction, with 
a Gaussian profile. 

We consider a 2D gas of N atoms rotating at fre- 
quency f2 around the z axis. The gas is confined in a 
harmonic potential, with frequencies ui x = wyl — v 2 , 
uj y = uiy/l + v 2 along the x,y axis respectively. The 
state of the gas is described by a macroscopic wave func- 
tion ip normalized to unity, which minimizes the Gross- 
Pitaevskii energy functional. In the following, we choose 
lu, huj, and y/h/JrntJ), as units of frequency, energy 
and length, respectively. The dimensionless coefficient 
G = Na s /a z characterizes the strength of atomic inter- 
actions (here a s is the atom scattering length and a z the 
extension of the wave function in the z direction for the 
initial 3-dimensional problem). The energy in the rotat- 
ing frame is 
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where Hn is defined by 
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and L 2 = i{yd x — is the angular momentum. We 
are going to study the fast rotation regime where fl 2 ap- 
proaches the critical velocity £l 2 := 1 — v 2 from below. 
Thus, we define the small parameter e by e 2 = \ — v 2 — Vl 2 . 
The spectrum of the Hamiltonian ^ has a Landau level 
structure. The lowest Landau level is defined as (see [l4j |) 

f[x + if3y) e [-^^ +{l3v)2 )\- l ^ x \ / is analytic (3) 

where 7 and (3 are some constants related to f2 and v 
given in the appendix; (3 is close to 1 if v is small. For 
such functions, < Haip,ip > can be simplified (see the 
appendix and [l4j). and in the small e limit (with e <C v) 1 
we are left with the study of 

ElllW = !\ (e 2 x 2 + K 2 y 2 ) |V| 2 + ^fdxdy (4) 



where n 2 ~ {v 2 + e 2 /2)(2 - ^ 2 )/(l 



This energy 



only depends on the modulus of ip. Hence, it is possible 
to forget the phase of ip, and use a simplified definition 
of the LLL: 

iP(x,y) = fix+i^e^wi^+^y 2 )}, f is analytic. (5) 

We recall that the orthogonal projection of 
L 2 (R 2 ) onto the LLL is explicit 0: U LLL (ip) = 

^ $ e -^ 2 - 2zTr+ \ z 'f)i>(x' ,y')dx'dy' , where 
z = x + iff y a nd z' = x' + «/3y'. We refer to the 
appendix of [16[ for details on the operator Ulll, its 
kernel and the computations: if an LLL function ip (i.e 
V> satisfies (JS])) is the ground state of ((H), it is a solution 
of the projected Gross-Pitaevskii equation: 



n 



LLL 



= 0, (6) 



where p is the chemical potential. 

The ground state of (0} without the analytic constraint 
is the inverted parabola 
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where i? T 



;^) 1/4 , ^ = (f^f) 1 ^ . Note that in the 
isotropic case v — (that is k = e), one recovers the 
standard circular shape R x = R y = [AG /(ire 2 )] 1 / 4 . Since 
k ^> e, R x is always large. On the other hand, the be- 
haviour of Ry depends on the respective values of e and 
k ~ v\[2. We find that R y is large if v <C e 1 ^ 3 while 
R y shrinks if v e 1 ^ 3 . We are going to see that in the 
first case, the profile ((7|) is reached in the fast rotation 
limit in the LLL using a vortex lattice, exactly as in the 
isotropic case, while in the second case, ([7]) is not a good 
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description of the condensate because the properties of 
the LLL prevent R y from shrinking, and in particular the 
energy is much higher than that of 0. 

In the first regime v <C e 1 / 3 , which we call the weakly 
anisotropic case, figure [T] provides a typical vortex config- 
uration, together with the corresponding density plot. It 
is obtained by minimizing the energy as a function of the 
location of vortices z\ with a conjugate gradient method. 
The vortex lattice can be described as in the Abrikosov 



FIG. 1: An example of (a): a configuration of the zeroes (b): 
density plot. There are 58 vortices with 23 visible vortices. 
v = 0.03, fl = 0.9985, e 2 = 2 x 10~ 3 , G = 3. 



problem 0, H| using the Theta function: 



[x,y,T) = 
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T/7 
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where z — x + i0y and t = tr + irj is the lattice pa- 
rameter. The zeroes of the function <fi lie on the lattice 

(Z © It) and \<p\ is periodic. The optimal lattice, 

that is the one minimizing b(r) = J \<j>\ 4 /(J \<fi\ 2 ) 2 is tri- 



angular, which corresponds to r 



a 2m/3 



(the integrals 



are taken on one period). As in the isotropic case [8[, 
we can construct an approximate ground state of ([4]) by 
multiplying the solution ([5]) of the Abrikosov problem by 
a profile p varying at the same scale as ptf defined in ((7| • 
Since this product is not in the LLL, we project it onto 
the LLL and define v = TLlll {p(%> y)4>( x i V\ T )) whose 
energy is 



p- 



Gb( T ) 



p 2 dxdy, 



up to an error of order y/ne^n 3 /e) 1 / 8 . Then, min- 
imizing with respect to p yields that p(x, y) = 

1 "'~ ( 6(7774' MVW4) where Ptf is given by © 



=Ptf 



The condensate indeed expands in both directions, and a 
coarse-grained density profile is close to the anisotropic 
inverted parabola. The vortex lattice is not distorted 
by the anisotropy since is close to 1; it is still trian- 
gular, as displayed in figure Q] Nevertheless, as in the 
isotropic case 0, H| , the lattice is distorted on the edges 
of the condensate, thereby allowing for a coarse-grained 
Thomas-Fermi profile in the LLL description. 
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FIG. 2: An example of (a): a configuration of the zeroes (b): 
density plot. There are only invisible vortices (32 vortices). 
Here, v = 0.73 = 0.6820, e 2 = 2 x 10" 3 , G = 3. The 
extension in the y direction is given by (|9} 



However, when v ^> e 1 / 3 , that is for fast rotation, the 
behaviour is very different as illustrated in Figure [2j we 
are going to see that the ground state is close to a Gaus- 
sian in the y direction multiplied by an inverted parabola 
in the x direction. There is no vortex lattice. There are 
only invisible vortices whose role is to create the profile 
in the LLL. The function Q does not provide the cor- 
rect behaviour of the ground state: though R y in is 
small, the condensate does not shrink in the y direction 
but keeps a fixed Gaussian profile [It} 



9(x,y) = — 
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We are going to prove that if u is the ground state 
and p(x) its projection onto the Gaussian ([9]), then p(x) 
is almost an inverted parabola, fndeed, in the LLL, 
we have the key identity J{d x \i)\) 2 + {1/ f3 2 ){d y \ip\) 2 = 
(7/4/3) J \ip\ 2 (see @). By adding and subtracting 
k 2 /(27/j) J \ip\ 2 to the energy, and using this identity, we 
find 
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This expression of the energy allows to analyze separately 
the contributions in the x and y directions. The ground 
state of ~{2/j 2 /3 2 )d 2 + (l/2)y 2 is the modulus of © and 
the ground energy is 1/(7/?). Projecting any function of 
the LLL onto the space generated by ([9]) times a func- 
tion of x, and using that e 1 / 3 jv is small, we find that 
E L ll(u) > (n 2 /{2 1 P)) + E 1D (p{x)) where 



E 1D (p) 




The minimizer of E\r> among all p's is of Thomas-Fermi 
type and we call it q: 



q(x) 



2 \ 1/2 




R 




since e 1 ^ 3 jv is small. This gives the energy estimate 

/ / \ 2 / 3 

mn,E LLL - — > E 1D (<1) - J- Q - ) . 

,( 12 ) 

Let us point out that this lower bound is optimal 
since we can construct a test function in the LLL with 
this energy. We project a Dirac delta function in the 
y direction times an inverted parabola in x, that is 
v(x, y) = An.LLL[5o{y)q{x)), where q is the function (TlTj) : 



v(x,y) = — e « u 



e *F 



^ x - x ' )2 - 2lx '^)q{x')dx' 



(13) 

The constant A = (2Tr/^(3) 1 ' 4 is a normalization factor. 
The fact that q varies on a scale of order e~ 2 ' 3 allows to 
expand (TIB")) in powers of e 2 / 3 : 



i—xy 



i-xy) (14) 



with an error of order e 4 / 3 . Inserting this expansion in 
the energy, we find 




le*a*q(x)* + ~^q(x)*j dx 

, (15) 
with an error of order e 4 ' 3 . This matches our lower bound 

Let us point out that according to (fTl"]) . the wave 
function v has no vortices in the bulk. This is corrob- 
orated by the numerical computation displayed in Fig- 
ure El Nevertheless, the inverted parabola profile in the 
x direction is obtained in the LLL thanks to the existence 
of invisible vortices, that is vortices outside the support 
of this parabola. 

Let us point out that the operator y 2 , whose ground 
state is the Gaussian Q is bounded below by a positive 
constant in the LLL: j R2 y 2 \ip(x, y)\ 2 dxdy > \ J R2 \ip\ 2 . 
This can be viewed as a kind of uncertainty principle [l8| . 
This decoupling in the x and y directions is possible only 
when the leading order term in the energy k 2 /(7/3) ~ v 2 , 
is larger than the energy of V Gve, that is when the 
ratio v 3 /e is large. When v 3 /s becomes of order 1, all 
the terms in the energy seem of the same order, the de- 
coupling in the x and y variables is no longer meaningful, 
and ([7]) does not provide the good behaviour either. The 
analysis in this intermediate regime is still open: it could 
display rows of vortices as obtained by [ll| . 
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The estimate of the energy (|12|) allows us to justify the 
validity of the model: indeed, the mean field approxima- 
tion is valid if the number N of particles is much larger 
than the number of one-particle states allowed by the 
chemical potential /z, that is N ^S> /Vmi- Thanks to (fl2|) 
and (I16|) . we find N > G 2 / 3 /e 1/3 

v. Since G is of order 
1, e 2 ~ 10~~ 3 and v < 10 _1 , this criterion is satisfied as 
long as N is greater than 10 4 , which corresponds to ac- 
tual values in experiments. However, if e gets too small, 
this condition gets violated and the states get correlated. 
The LLL approximation is valid if the ID energy E\d is 
much smaller than the gap /i 2 between the LLL and the 
first excited state: (Ge) 2 ^ <C 1. 

Conclusion: When the anisotropy is small compared 
to how close the rotational velocity is to the critical ve- 
locity, that is e 1 / 3 3> v, the behaviour is similar to the 
isotropic case with a triangular vortex lattice. A strik- 
ing new feature is the non-existence of visible vortices for 
the ground state of the energy in the fast rotation regime, 
that is when e 1 / 3 <C v. The profile of the ground state 
is a large inverted parabola in the deconfined direction 
and a fixed Gaussian in the other direction. Our analy- 
sis indicates that an asymmetric rotating condensate un- 
dergoes a similar transition as a condensate placed in a 
quadratic+quartic trap where at large rotation the bulk 
of the condensate does not display vortices (l9j]. Our in- 
vestigation opens new prospects for the experiments: in 
particular, if a condensate at rest is set to sufficiently 
large rotation, then vortices should not be nucleated. 



Appendix 

As computed in [3] on the basis of ideas of 
Valatin [20], the eigenvalues of the Hamilto- 
nian Hq are 1, yLi|, where /i 2 = 1 + £1 2 — 

x/i/ 4 + 4f2 2 , i4 = 1 + fl 2 + V^ 4 + 4ft 2 . We de- 
fine a = Vt/ 4 + 4f2 2 , Pi = (2fi/ii)/(Q!-2fi 2 + ^ 2 ), 
(3 = p 2 = (2ft/i 2 )/(a + 2ft 2 + ^ 2 ), 7 = (2a)/n, 



d = 
where 

i 

V2 



02 

-i\2dy - 



^ {-l\2d^dy 



( 7 AiA 2 )/2, c = (A 2 + A 2 )/2AxA 2 . Then 

= J f a l a l + a l a l) + 5 u4°2 + a 2 c 

^(-iAid-^x+cAiy) 
((iA] -1 — A 2 «i) x), and «i 
+ c\ 2 x) + ^= ((Aicd — d\ 2 1 )y — i\\d.. 



We have: 



»2, a 2 



= M2, 



ai,a{ 



Mi 



and all other 



commutators vanish. The LLL is defined by a 2 ip — 0, 

that is f(x + ip 2 y)ei-^(^ x2+ ^ (f32v)2 )h l £ xv , 
with / analytic. It is always possible to change the 
analytic function /(£) into /(£) ex P( — ^£ 2 ) m t ne above 
definition, since exp(— <5£ 2 ) is an analytic function of 
£. Hence, for S = v 2 / (8n/3 2 ), we find the alternative 
definition of the LLL ([3|), with /3 = fii- This defini- 
tion is equivalent to the one given by Fetter in [lil ]. 
However, contrary to [l4|], the coefficients in (J3j) are 
not singular in the limit e — ► 0. Indeed, in this limit, 
2 ~ ^/{l-v 2 )/{l-v 2 /2) and 7 ~ (4 - 2^ 2 )/ x /T _ ^. 
This is due to the addition of the above-mentionned com- 
plex Gaussian in the definition of the LLL. In the LLL, we 

have (H n ijj,ip) = \ (( a\a x + aia\\ ip, V>) + f V>) • 
We then express x and y as linear combinations of 



a\i a 2 , a{, & 2 



(HEl and get, if V € 



ViPix + jV z ) M'dajdj/ (16) 



A 2 = (a-2ft 2 + ^ 2 )/(2a), A 2 



(a + 2Q 2 + v 2 )/2a, 



which provides (|4|) with k 2 — 7 /ji/2/3i since Jfiipi ~ 2e 2 . 
Note that /ii ~ and /i 2 ~ (2 — ^ 2 ). 
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